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Abstract 


A  method  already  developed  for  the  diffraction  of  scalar  waves  by 
a  circular  disc  is  generalised  so  as  to  be  applicable  to  the  diffraction  of 
an  electromagnetic  wave.   The  generalisation  leads ^  in  a  straightforward 
manner,  to  an  asymptotic  development  of  the  current  density  on  the  disc  at 
high  frequencies.   The  principal  complication  which  arises  is  in  solving 
three  simultaneous  linear  equations  for  three  constants  which  occur  in  the 
method  of  solution.   Detailed  calculations  are  given  for  a  normally  incident 
plane  wave. 

Introduction 

The  problem  of  the  diffraction  of  an  electromagnetic  wave  by  a 
circular  disc  has  been  the  subject  of  much  investigation.   A  solution  was 
published  by  Moglich  [l]  as  early  as  1927.   However,  the  solution  which 
satisfies  the  edge  conditions  was  not  obtained  until  15  years  ago.   This 
solution,  which  is  expressed  in  terms  of  spheroidal  functions,  has  been 
derived  by  Meixner  [2],  Meixner  and  Andrejewski  [3],  Flammer  [h] ,   and 
Belkina  [5]  when  the  primary  field  is  a  plane  wave.   The  corresponding 
solution  for  a  dipole  field  has  been  found  by  Meixner  [6],  Flammer  [7] 
and  Belkina  [5].   Other  exact  analyses  have  been  given  by  Nomura  and  Katsiira 
[8]  and  Lur'e  [9].   Numerical  results  up  to  ka  =  10,  k  being  the  wave  number 
and  a  the  radius,  have  been  obtained  by  Andrejewski  [10,11]. 

Several  approximate  methods  for  tackling  the  problem  have  been  devised, 
some  are  appropriate  to  low  frequencies  and  some  to  high  frequencies.   At 
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high  frequencies  the  series  of  spheroidal  functions  become  unmanageable  so 
that  other  approaches  have  to  be  used.   One  method,  based  on  physical 
reasoning,  is  due  to  Clemmow  [l2]  and  Millar  [13,1^]  -  this  gives  the 
dominant  part  of  the  field.   The  method  of  Seshradi  and  Wu  [15]  is  based  on 
an  idea  of  Levine  and  Wu  [l6],  which  chooses  the  dominant  part  of  the  kernel 
of  an  integral  equation  so  that  it  can  be  solved  by  the  Wiener-Hopf  technique. 

The  object  of  the  following  paper  is  to  indicate  a  systematic 
process  which  provides  as  many  terms  as  one  is  prepared  to  calculate  of  the 
asymptotic  development  at  high  frequencies.   It  is  based  on  a  theory  which 
has  already  been  successfully  applied  to  the  diffraction  of  scalar  waves 
(Jones  [17] )•   The  generalisation  required  is  fairly  straightforward  because 
the  vector  integral  equation  of  electromagnetic  waves  can  be  reduced  to 
relatively  simple  scalar  integral  equations  (Jones  [I8]). 

In  Section  1  the  salient  points  of  this  reduction  are  repeated. 
Section  2  is  concerned  with  converting  a  typical  integral  equation  into  a 
singular  integral  equation  of  the  first  kind  which  is  especially  suitable 
for  handling  at  high  frequencies.   From  that  point  onwards  the  analysis  is 
specialised  to  a  normally  incident  plane  wave.   Only  two  integral  equations 
have  then  to  be  solved.   The  asymptotic  solution  of  one  is  derived  in  Section 
3  and  of  the  other  in  Section  k.      The  solutions,  so  obtained,  involve  three 
constants  which  have  to  be  chosen  so  that  the  field  behaves  in  a  certain 
way  near  the  edge.   These  constants  are  found  in  Section  5-   A  general 
formula  for  the  scattering  coefficient  is  derived  in  Section  6  and  its 
evaluation  for  high  frequencies  is  given  in  Section  7- 


In  order  not  to  interrupt  the  argiiment  of  the  main  text^  various 
calculations  have  "been  carried  out  in  appendices.   References  to  equations 
in  appendices  are  written,  for  example,  (Al),  meaning  equation  (l)  of 
Appendix  A. 

1.  The  integral  equation 

Let  all  lengths  be  normalised  so  that  the  radius  of  the  disc  is  unity 

2    2 
and  choose  coordinates  so  that  the  disc  occupies  x  +  y  <  1,  z  =  0.   Let 

the  incident  plane  wave  have  electric  vector  E  with  Cartesian  components 

given  by 

-ia(x  sin  9  +  s  cos  9  ) 

E   =   (cos  CD  cos  9  ,  sin  cp  ,  -  cos  cp  sin  9  )e  , 

-o  o      o      o       o      o 

the  time  dependence  e    being  understood  and  omitted  throughout.   The  plane 

wave  is  propagating  in  a  direction  whose  polar  angles  are  9  /p  ;  in  the  case 

of  normal  incidence  9=0  and  cp  =  0.   The  constant  a  is  the  product  of  the 

o  o 

wave  number  k  and  radius   a  of  the  disc;   it  will  be  taken  as   a  positive  real 
nianber  throijighout. 

The  total  electric  field  at  the  point  R  can  be  written  in  the  form 

(1)       E(R)   =  E^(R)  +  (grad  div^  +o?)   j    l{\)   t(iR,%)dS 

S 

where  S  is  the  unit  circle  on  z  =  0,  R  is  a  point  of  S,  t(RjRi  )  =  e    ' i'/IR'^^U 

div  is  the  two-dimensional  operator  independent  of  z  and  J  is  a  constant 
multiple  of  the  current  density  induced  on  the  disc.    Since  the  disc  is 
perfectly  conducting  the  tangential  components  of  E  vanish  on  the  disc  and 


so  the  integral  equation 

(2)       E^^  +  (grad^div^  +  cx^ )    JiE-^HiE,E-^)^^     =  0     (R  in  S] 


is  obtained  to  determine  J. 

The  integral  equation  (2)  can  be  converted  into  a  set  of  integral 
equations  for  the  Fourier  components  of  the  current  and  charge  densities 
(Jones  [l8]).   Recently,  other  approaches  have  been  given  by  Lur'e  [9] 
and  Williams  [19].   Let  p  =  div  J  and  let  r,cp  be  polar  coordinates  on  the 
disc.   Since  p  and  J  are  periodic,  with  period  2-n,    in  cp  we  can  write 

(3)  oo  0°  c° 

n=-oo  n=-oo  n=-oo 

Then  it  has  been  shown  (Jones  [l8])  that  p  ,  a  and  b  satisfy,  for  0  <  r  <  1, 

(h)  1 

1     .  p 

-i   tan  G  cos  9  e'^'^^J   (a  sin  9^)  +  A  J   (ar )  +  a   /   p    (r    )t   (r,r    )r  dr        =  0, 
o  o  n  onn  'nxn±x± 


o 


(5) 

2rt 


i;    /^^or-'-^^ocp)^"'''^^^-^!:-?)   /pn^^l^^n^^'-l^-l^-l^^'/^n^-l^Vl^^^-l^^l^^l  =  '' 


(6) 

2n  1                                                        1 

t    /(Eor-^^ocp^^"'"'^<^^n^^S)  /Pn(-l)^n(^'^lK^-l^^'/'n^^l^Vl^^'-l^-l^^l=° 

o                                              ^        '^  o                                                     o 


where  A  is  an  arbitrary  constant  to  be  determined  eventually  by  the  edge 
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conditions  and 

2        2  2 

2jr     -ia(r  +  r^   -  2rrT  cos  9) 

^— T cos  n9  dcp. 

(r    +  r,  -  2rr, cos  9) 
oil 


Substitution  for  p  from  (k)   in  (5)  and  (6)  leads  to  the  integral  equations 


(Q)     e'^^^'^'^'^^icos   9  sec  G  +i  sin  9  )J  ,  (a  sin  9 


o'  n+1' 

1 

+  A^J^^^(ar)  +a2j^a^(r^)t^^^(r,r^)r^dr3_  =  0, 
o 

(9)  e"^^"~"'"^'^^(cos  9  sec  0  -i  sin  9  )J^  -,  (a   sin  Q    ) 

'     o     o        o  n-±       o 

1 

-  A  J  -,  (ar)  +a    /  b  (r-,  )t  n  (r.r-,  )r-,dr-,  =  0, 
n  n-1^  '  J        n^   1'   n-1^  •*  1  ^  1  1    ' 

o 

when  the  specified  form  of  E  is  used  in  (5)  and  (6).   In  the  particular  case 

of  normal  incidence  (0  =  0,  9  =0)  the  incident  wave  term  disappears  from  (h), 

is  absent  from  (8)  unless  n  =  -1  and  from  (9)  unless  n  =  1.   Since  there  will 

be  no  field  if  there  is  no  incident  wave  we  see  that  a  =  0,  b  s  0  (n  ^  ±l)  and 

n       n    •    ' 

that 

(10)  1  +  A_^J^(ar)  +  a^   /  a_^(r^)i|r^(r,r^  )r^dr^   =  0, 

o 

1 

(11)1  -  A^J^(ar)  +  a^   r  b^(r^  )tQ(r,r^  )r^dr^  =  0. 


Since  A_^  and  A  may  not  be  zero  we  must  also  consider  the  equations  for  a^  and 
b_^  namely 
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1 
(12)      A^JgCar)  +  a^  /  a^(r^)t2(r,r^)r^dr^  =  0, 

o 


(13)     -A_^J2(ar)  +  a^  /  b_^(r^  )t2(r,r^  )r^dr^   =  0  , 


since  i|^  p  ~  ^o*   Thus,  only  p  and  p_  will  be  non-zero;  a  fact  which  is  also 
evident  from 

P   =     -:r     T~    (^(a  +  b  )}  +  :7-(a  -  b  ). 
n     2r  dr     n    n     2r  n    n 

At  the  edge  the  radial  component  of  the  current  must  vanish  (Jones  [20])  and 

this  requires  a^  +  b  =  0,  a_   +  b_   =  0  at  r  =  1.   From  (10)- (13)  it  is  immediately 

clear  that  A  =  -A   and  that  a   =  b  ,  a   =  b   .   Hence 

1   d  ^^1 

•^r  =  ^^l"'^l^*'°^  """^   J^  =  (a^-b^)sin  cp   and   p=-  -^   i^^  +  \)   +  —   • 

The  further  edge  conditions  that  J  and  p  be  integrable  will  be  met  if  neither 
a^  nor  b,  is  more  singular  than  (l-r)  ^  as  r  -^1. 


2.  The  transformed  integral  eqiiation. 

All  the  integral  equations  {h) ,    (8)- (13)  are  of  the  same  type,  namely 

1 
(lU)  /  r^k(r^)t^(r,r^)dr^  =  g(r)  (0  <  r  <  1) 

o 

where  g(r)  is  of  specified  form  and  k(r,  )  is  to  be  found.   To  convert  this  integral 

equation  to  one  of  a  more  convenient  type  we  use  the  result,  proved  in  Appendix  A, 

(A4) 


n 


*n(r,r^)  =  2rt(- 


.n       r 


n-1  n 
a   r 


.rdr  y 


miiHr,r  , 


(r  -  X  )  ' 


P   pi       P   p  i   (^1  -  ^' 


2   2,i(n4; 


i  /  Vi^°'^''l"  "^  )'^Ji^^(^  -^  ^"^  7^ 2^711 


dx 


r  -  X 


Thus  (1^)  can  be  written  in  the  form 


(15)   2rt(-) 


n  r 


n-1  Vrdr 


a 


r  T   .  ,2       2A, 


■-  o 


,~2        2.1/1+ 
(r  -  X  J  ' 


h(x)dx 


p    Jila(r  -  X  )  , 

'  ~f2        2^1/1| 
(r  -  X  j  ' 


hfx)dx 


=  g(r; 


in  0  <  r  <  1  where 


(16)  h(x)  = 


.    ,   2   2.^.,   2   2.i(n-i),  .   s  1-n^ 
J  iia(r,  -X  )]{t^-x    Y^      '^'k(r  )r,   dr.,. 
n-2    -L        X  J.  X    X 


X 


There  are  two  methods  of  tackling  (15).   One,  which  leads  to  equations 
similar  to  those  obtained  by  the  author  earlier  (Jones  [l?]  ) ,   will  be  given 
here.   The  other,  which  leads  to  equations  analogous  to  those  of  Bazer  and 
Hochstadt  [21],  is  described  in  Appendix  B. 


Let 


1  f   2   2.n-l 

(17)         ^(r  )   =   /   — -7  -rri  ^^ 

^      -*   (n-l).'2'^  ^     u^  ^ 
r 
1 


for  n  >  1  and  k  (r  )  =  k(r  )  for  n=  0.   Then 


■    .   1    1  ^  >/  2   2.n-2 
^--J       ^^^^^   =  fe^    ((n-2):2-^    .-^" 


1 
__i._   r  (_)n-l  k{ul  ^^ 

^1 


(18: 


:-f^{r^)/rl   . 


Note  that  because  k(u)  is  no  more  singular  than  (1-  u  )  ^  as  u  -> 1,  (1?) 
shows  that  K(r  )  is  proportional  to  (1-  r  )  "^  as  r  -^1.   Further,  since 
k(u)  is  bounded  at  u  =  0,  K(r_  )  is  bounded  at  r,  =  0. 

On  account  of  (l8),  (l6)  can  be  written 

1  ,    .  n 

X 

For  n  =  0  this  formula  need  not  be  modified.   For  n  >  1,  integrate  by  parts 

in  such  a  way  as  to  reduce  the  derivative  on  k  .   Because  ^  77  '■^  J  (z)}=  z   "^v_i( 

we  have 
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h(x) 


J  i.la(r,  -  X  )  j  (r  -  X  )  ' 
n-2    -L         -L 


;-r 


n-l 


-il 


^(r^) 


X  ^^"2 


"  1  ly 


=  a 


r  X    r  ^  2   2.iw  2   2.i(n-j 
X  '^'  2 


^(-'""'(dr^"'»(-i)^i 


since  the  first  term  vanishes  because  of  the  behavior  of  the  Bessel  function 
near  r  =  x  and  of  k  near  r-^  =  1.   The  remaining  term  is  a  constant  multiple 
of  the  original  integral  but  with  n-l  for  n.   Repetition  of  the  process 
leads  to 


(19) 


h(x)   =  a 


n  J  i{a(r  -  x  )^) 

/  _=2 r  K(r^)dr^ 

^      /  2   2s7 
x     (r  -  X  )4 


This  is  also  valid  for  n  =  0  in  view  of  the  remark  above. 

Substitute  (19)  in  (15)  and  interchange  the  order  of  integration  in  the 
second  term  on  the  left-hand  side.   The  result  is 


2.(-A^(4j 


J_i.{a(r^-  x^)2}  -1  T  ,r./^2  ^2,k 


(20) 


{    (r^-x^)^/^  J 
o  ^  X 


2     2^3■ 


— ^ r, r^  K(r-,  )dr^ 

T2   2.l/i|     1  '  1'   1 
(r-^-  X  )  ' 


2   2. 


J  i{a(r^-  x^)2}  Ji(a(r^-  x^)^) 

_  2      -L  P 


'/"l^^^l^/     /2   2.lA     ^2   2,1 /U   ^"'^l 


.2  iTi/T 

(r^-  X  )  ' 


,~2       2,lA 
(r  -  X  )  ' 


=  g(r). 
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But  it  has  been  shown  elsewhere  (Jones  [17])  that 


p^   J_l(a(r^-  x^)^]   J^(a(r^-x^)^3         o  J_i{a(r^-  x^)^}   Ji[a(r^-  x^)^] 
^     (r  -  X  )  '       (r  -  X  )  '  ^     (r  -  X  )  '      (r  -  x  )  ' 


Therefore  (20)  can  he  written  as 


2   2^in 


/   -vn   p  J_l{a'(r  -  x  )^} 
(21)      2,(-)n,n^Q^j   J  -^^2_^2^lA   ^^^^^^   =  ^^^^ 


where 


o 


^  J  i{a(r^-x2)2}  /^  Ji{a(r^-  x^)^: 


(22)  X(x)   =   /    ^  2   2,lA'  r^-(r^)dr3_-i  /  -^^-^3711-  ^l^^^l  ^^1 
'^     (r^  -  X  )  '  ^         (r^  -  X    )   ' 

X      1  o     1 


Integration  of  (21)  n  times  gives 

f  '-'^'f-  fp  ,(,),,   .  y.r^.-J rV  ,2,„-i  at^  ,„ 

J     ,  2   2vl/l|  /  I  m      ,   ^,,^n    J  ^  n-1 

'-'     (r  -  X  )  '  ^ — r        (n-1). 2  rex  ^  u 

o    ^  m=0         ^  o 


where  c  ,c  ,...,c  _  are  arbitrary  constants.   If  n  =  0  the  right-hand  side 
is  replaced  by  g(r)/2rex. 

Replace  the  Bessel  function  by  its  trigonometric  equivalent  and  then 

r  2   2  -            n-1                       r 

r  cosfO'^r  -  X  )^1  ,,  ^^      'X"'  '  2m   1 P  ,   2   2. n-1  gfu)  ^ 

/  :r' o  1  —  X{x)dx     =   >   c  r   -  _,!  '1  1  /  (u  -  r  ^~~7   du 

J  /  2   2.2              /  I  m  ,      ^  Mo^+s  2  2  J  n-1 

"^^  (r  -  X  )               ^— 7'  (n-1)  2   rt%^  ^  u 

o  ^                     m=0  ^  o 

1  1 


^/2  ~  - 
If  n  =  0  the  right  hand  side  is  g(r)/2  '  jt'^a^ 
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The  solution  of  this  integral  eqiiation  is  known  (Jones  [22] 

2 


to  be 


(23)   A(x) 


1  _d_  f     cosh(a(x  -  y  )^  ] 
rt  dx  J     ,2        >2 


m 


n-1 

L  m=0 


(n-l):2  ^jr^a 


'^   t^t  o 


y        1 

/  (u-y)      ^;r  ^u 


dy 


If  n  =  0  the  term  in  [   ]  is  replaced  by  g(y"')/2  jt^a^. 

Row  that  X(x)  is  known  K(r  )  is  determined  from  (22),  which  may  be 
expressed  as 


cos  (a(r  -  X  )^  ] 

~~2       FJ 
(r  -  X  ) 


p    sinia(r  -  x  )  J 
r^K(r^)dr^  -  ij  ~--       ^ 


r^fc(r^)dr^  =  (-na)^X(x; 


2       2 
Take  as  new  variables  z  =  x  ,  t  =  r  and  the  equation  becomes 


(24) 


1  1  i 

^°^^^(^:^f)    ,(ti)dt  -  i  r  ^^Mt^  ^^^i^^^  ^  (2::a)2x(z2; 


1    1 


(t-z)' 


t-z 


valid  for  0  <  z  <  1.   In  this  foimi  the  integral  equation  is  of  the  same  type 
as  (6)  in  our  earlier  paper  (Jones  [l?])-   Hence  we  deduce  that 


(25)   /  wk(w; 
o 


-ia(w+v)    -ia(w-v) 


w  +  V 


w  -  V 


dw 


K^°^(v) 


(0  <  V  <  1) 


where 


(26)   K^°^(v; 


2 

V      1 


/I   >'2  d 


-  [z: 


dv 


X(z^)cos(afv^-z)^l 
(v^-z)^ 


dz 
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and  the  principal  value  of  the  integral  in  (25)  is  to  be  taken.   Subsequently, 

a  principal  value  will  always  be  assumed  where  necessary. 

A  solution  of  (25)  is  required  that  has  no  more  singularities  than 

the  solution  of  the  corresponding  integral  equation  in  our  earlier  paper 

so  that  the  same  method  as  was  used  there  may  be  adopted  here  to  find  an 

asymptotic  development  for  large  a.   This  asymptotic  development  will  involve 

the  constants  c'  because  of  (23)  and  these  constants  will  have  to  be  chosen 
m  ^ 

2  n-- 
so  that  K;(r,  )  behaves  like  (1  -  r  )   ^  as  r  ^  1.   Once  k   is  known  k(r  ) 

can  be  determined  from  (l8).   Thus  the  distribution  on  the  disc  will  be 

known  and  all  other  quantities  can  be  deduced  therefrom. 

In  order  to  illustrate  how  the  details  work  out  in  practice  we  shall 
now  consider  the  particular  case  of  the  plane  wave  at  normal  incidence. 

3.  The  plane  wave  at  normal  incidence. 


It  has  been  pointed  out  in  Section  1  that,  for  the  plane  wave  at  normal 
incidence,  it  is  necessary  to  determine  only  a,  and  b  .   In  this  section  we 
consider  the  determination  of  b  .   Equation  (11)  is  of  the  same  form  as  (ll|) 
with  n  =  0,  k  =  b  and  g(r)  =  (A  J  (ar )  -  l}/a  .   Hence  the  theory  of  the 


-ia(w+v)    -ia(w-v)  /  -. 

wb,  (w)<^^ +- V  dw  =  B^°^v, 

1      1    W  +  V  W  -  V 


B^°\v)  =  (, 


1   vi  d   r  cos(a(v  -z)^]  2  i  d    fcoshfaf  z-y)^  ] 

'^   ^   dv  ,/    ,  2   v2      3t    dz   /     ,    ^2 
"-'     f  V  -7.  )  ^      1  z-v 


(v  -z) 


(z-y)^ 


A^J^(ay^)  -  1 
23/2  A5/2 


dy  dz. 
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Now 

i  z  _ 

,2  ■ 


r_    /     ^^   cosh(a(z-y)^}    ,  lim  T  t   /        2n    cosh(a(z-y)'  j 

,/o  /         ^2  a— >a/0-L  /         ^2 

"^c.  (z-y)  1  ^Q  (z-y) 

1  ^^  (z-y) 

211^  3 

■^™  /   J   (a-,z^sin  0)(2jta)^I    i  (az^cos  9  )z   sinG  cos^GdQ 

a     -*  a    J      o  ^    1  ^  -2 

o 

2 

on  putting  y  =  z   sin  0.      This  last  integral  is  a  particular   case   of   Sonine's 

second  finite  integral    (Watson    [24] )   and  so 


o 


1  i.r/^2-1  n-  i|ii    Ji{z^(a^-a    )^} 

T   /      2\   coshia(z-y      j ,  lim  ,^      ^2    4,      2\-2     ?  1 

Jjay    )  -    /    ^    ■'  '     'dy  =  ^^  (2na)2z    (az^ )        -^j ^— 

(z-y)  1  i;    2     2,1; 


z    (a   -a    ) 


•        i^    2     2.i 
sm   z    (a^  -a    ) 
lim         p  1        ' 

a      ^  a  ,2  2^2 

1  (o^-L  -  a.    )'' 


(28)  =        2z 


Hence 


2 

2_N2T  i 


^(o),    V               1       d        f    cos(a(v  -z)^l    .„  ^       i,^ 

B'      (v)      = ^—      /     ^ — i.    '    '    (A^   -   cosh  az'^)dz 


2jn.     "'     ^  (v^-z^^ 

o  ^ 


^    {4a.  cos  av   -    4  cos  av   -    4xv  Si(av)) 

4jta 


-  Ill  - 

X 

where  Si(x)  =   /  — 7 — dt. 
o 

The  last  two  terms  in  the  {   )  of  B    are  the  same,  apart  from  sign, 

as  those  which  arise  for  the  scattering  of  a  scalar  plane  wave  by  a  sound-soft 

disc.   The  part  of  b  due  to  these  terms  may  therefore  be  deduced  from  our 

earlier  paper.   Denote  by  f(w)  the  solution  of  (27)  when  B^   (v)  =  4  cos  av 


+  k  cLv   Si(av).   Then 


(29)  b^(w)  =  -—  {A^d(w)  -  f(w)} 

kztCL 


where 


'  -ia(w+v)    -ia(w-v' 

(30)      /  wd(w)  { +  }     dw     =      k   cos  av. 

^       '  '         ^     '   ^      V  +  V  w-v 


O  ^  _y 


Let 

/-n  \  T^/  \      r  wdfw)   -ia(w-v), 

(31)  D(v)  =   /  — ^ — ^  e       'dw  ; 

^       J   w  -  V 

o 

then 

1      1 
law     p   -,    2      -r,/    \    -lav 

(32)  d(w)   =  -f^ r   r  (i^)    ^^^ dv 

A^(i-w)^  J      ""         "  -  ^ 

^       o 


1       .  ^ 
1   -lav   p  ^,,   s  -lat   -,  ,  1 

^L.f  e /  D^tie fi^)2dt 

1+v  ^    jt    ,  /    t  +  V     t  ' 


as  in  our  earlier  paper.   The  eq^uation  satisfied  by  D(v)  is 
D(v)  =  k   cos  av  -  ( 

_1_  I  V  Jl 

o 

for  V  in  (0,1).   Write 

1 
D(v)  =  2e    +  (^)  e    G(v) 

and  then 
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»M --i/fey-^f¥^- 


This  equation  is  analogous  to  equation  (31)  of  our  earlier  paper. 
There  is  the  difference  that  2  replaces  h, (v).   However,  the  argimient 
which  followed  (31)  can  still  be  carried  through.   Hence  we  can  form  the 
iteration  scheme 

n^  ^       V    ,/    t-v 


00  1 

lay  p  2 

V^)   =  ^^  /  ^ifc)    e^(x)M(x,y)dx       (n>l), 


o 


^t-i*  Vi^^)^ 

1 


-2iat 


e>)   =     (^^-)  -^^^^i-rr^ dt  (n>2), 


(33)      e^(v)   =  2 


and 


i^xyHJ^^(ax)J  (ay)  -  Ify^H?^^ ''(ax)J,  (ay ) 
M(x,y)   =  i °  o       1 


2    2 
X   -  y 


From  the  general  theory  of  our  earlier  paper  we  can  say  that  each 
iteration  is  smaller  than  the  preceding  by  the  factor  a  ,  as  far  as  a  is 
concerned.   Also,  from  (32),  the  asymptotic  development  of  d(w)  is  given  by 
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,/    N  e /  /1-v^  2dv 

d(w)  ^  -  — r— [ r      /  ( )  

2    2/-,         \2     J  V  v-w 

^  o 


^    /    V    -law  .  00  1    -c    / .  -i    -2iat 

E    ( w  )  e  law  n    ,    ^    2  E    ( t  j  e 

-^- -.-4 r      fi^)     -^, dt 


(31,)  n=l  I  ^   (1-^^        1 


rvj  r 


2^iaw 


^  ^    .     >    -law  .  00  I   rr    f+\    -2iat 

E    ( w )  e  law  n    ,    .,    -?   E    ( t  j  e 

r  +   >  l-"^ +-? r  f{^)     ^ dt 


n=l 


The  sulDstitution  of  (B^j-)  in  (29)  together  with  the  known  form  of  f  now 
gives  b  completely  apart  from  the  unknown  constant  A.^  . 

k.    The  determination  of  a  . 

The  integral  equation  for  a,  is  (12).   This  i?  the  same  as  (ik)   with 

the  identification  n  =  2,  k(r  )  =  ^1(^1  )  ^^'^  si^)   =  -A,Jp(ar)/a  .   According 

to  the  theory  of  Section  2,   and  (I7)  in  particular,  we  introduce  the  function 

p  (r  )  defined  by 

12   2 
o  u  -  r 

(35)  Pi(r^)  =  2  J  ~^I a^(u)du 

^1 

so  that,  on  account  of  (I8), 


(36)  a^(r^)   =  rl     (;^)\(r,: 
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1 
2.2 


The  function  p  (r  )  must  be  bounded  at  the  origin  and  behave  like  (1-  r  )  as 


r.,  ->  1. 


The  function  corresponding  to  X{x)   is,  from  (23), 
2 


X         2   - 

1  d   r  coshlafx  -y)^3 


jt  dx  ,/      ,2      ,t 

X.        (x  -y) 


A. 


c'  +  c'y  + 


1 


-o   ^l-'    I  ^ 


-^  J^(au2)du 
u    2 


o2  2  t  o 
2  a  Jt 


Now 


(— )   Jjau)  =  cc  —^ 


u 


so  that 


y 


-^  J  (au2)du 
u   2 


=  2 


■u  -y 


u'  '2 


)j  (au)du 


^  Ju{u^-y)    (^-)  Jjau)du 

a   "-^ 
o 


a 


(u  -y)  — —  J  (au) 
^    -^  ^  udu  o^   ^ 


1  i 

2  2 

ly  y 

-  ^  77.   J  (a-u)du 


2  ,/   du  o' 


a 


-y  +  -|  (1  -  Jo(ay2)) 

a 


Hence  the  function  corresponding  to  X{x)   is 
2 


1  d    r   cosh(afx  -y' 


rt  dx 


2   ^t 


(x  -y)' 


d^  +  d^(y  -  2/a  ) 


A^J^(ay2) 


(Scc^jt 


9.^t 


dy 
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where 


c'    +  2d, /a     +A,/(8a^jt)2  H      =  c'    -A,/(l28a 

o  1'  1'  ^  '    '  1  1  1'  ^ 


\9 


On  account   of   (28)  the   fimction  corresponding  to   >i,(x)   is 


1  J               ^'^l                           ^1                                 \ 
—  ■\2(d    -  — -)   cosh  ax  +  x   sinh  ax — 


rt   1     '    o  2 

a 


Therefore,    from   (25)  and   (26), 


1 


ia(w+v)  -ia(w-v)l  ,    ^ 


;37)    j  -Pi(^M'^~"T-  ^  ^— 

o 


where 


Pi°^(v) 


/^^2    ^ 


2        i 
^2  jiL.      [    COS    {a(v    -z)^] 

2jt^      dv      /  "2         ^i 

^  (v   -    z  ) 

O  ^ 


2d^  i.        ^-L     i  i  \ 

X    •(2(d    -  — ;::)   cosh  az^   +  z^sinh  az^    -  FTT?'     ^2, 

°       a^  ^  (2m9)2l 


Now 


"1 


Pr'(v)      =      (^)^(d  1 


2d^  i  2A  cos  av 

,    ,^  )    {\  cos  av   +   iptv   Si(av)}    -    {f-^  ^^ 

a  (2na-^j2 


^2jt^     dv 


8d  "2 


1^~  i  2       i 

— r-  sinh  az'^sin{a(v   -z  )^  } 


TV 


a 


1 


llrl  i 

/a  v2     d.        r    1      .     ,    ,   2.      >2-,   J  sinh  az^    ,  ,         2 '^    ■, 

('5")     ^      /      — ^  sinla(v   -z)    J  < T +  a  cosh  az    >  dz 

o 
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a 


,4  cos  av   +  kav   Si(av)} 


^  2A  cos  av 


(2rta^)' 


^2rt^ 


^^v  f    ,    2      .^^ 
1      cos  ia(v   -z  J    J 


a 


1. 

sinh  az^    ,  ,  i    •, 

— I +  a  cosh  az^      dz 


V 


2d  2A  cos  av        8d  v 

'd    -  — ;:;)    [h  cos  av   +   ijav   Si(av))    -  r— t"  + Si(av 

°        a^  (2Tta^)2 


kd 


1  2   2 

+  {sin  av   +  av   cos  av   +  a  v   Si(av)) 

a 


^2-a' 


2A  c  OS  av 

d    [k  cos  av   +   ijav   Si(av)} ,-,   i 

(2m9)2 


1  2   2  '^3 

+     — -      [-2cos  av   +  av   sin  av   +  a  v^cos  av   +  a'v   Si (av ) 

a 


Consequently,  remembering  (30),  we  see  that 


(38)  p^(w: 


{^J^  <|d^f(w)  -  -^-TTX   +  d^p(w; 


'■2-e'      )  o 


9^2 


(Bm^) 


where 


wp(w) 


-ia(w+v)      -ia(w-v) 


w  +  V 


w  -  V 


dw  =  P^°Vv 


and 
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(39)      P        (v)    =  -^    (-2   cos  av  +av   sin  av+  a  v   cos  av  +a  v   Si(av)], 
a 


Let 


■wp(w) 


-ia(w-v ' 


w   -    V 


=      P  V. 


Since  p  ,  f  and  d  are  all  boionded  at  the  origin  so  is  p. 

1 


Hence 


law 


(1+0)    p(w)  =  -  -^^ T 

«  W2(l-w)2 


V        V  -  w 


and  the  equation  satisfied  by  P(v)  is 


P(v)  =  F^°\v) 


i   -iav 
^l+v'^      It 


^  t  t  +  V 


Put 

00         ,      , 

p   -ia(t+v) 
P(v)  =  2e^'^''(aV-iav-  2)^+  2iav3   /   ^  ^  ^  ^  '  dt  +  P^(v) 

-2v 


and  use  the  results 


1  1 


2,2 


1  i 

2    r /i-tx2  ;  2, 


(1^)"   2(a^t   -iat-2;   at   =  ^   /    {^f  ^a^t   -  a(i+av)    +  ^'-^ 


2   2 
a  V   +  iav  -  21 


dt 


a    (t   +  v! 


a 


(k-l) 


2n 


i  a^   -  J  a(i+av)+  (aV+iav-2)    [{^f-   1} 


(lzW)2   2ia^    /  e ^^  ^^  ^ 


■  v    J    t   -  w 
o 


-2ia    /   e 


1  1 

-iat    r /1-w  2 


(~)  r^^""  V+t)(w+v)'^(w-t)(t+v) 


21 


-2ia        e 


-iat 


8  '^  +i«(^-^)  -^TT^V"^ 


-  1>  + 


t\ 


t+v      '    t 


^)"H(t-l)-   1 


dt 


from   (ElO)   of  Appendix  E  of  Jones    [l?] .      The  last  integral   can  be  written 

P  j^    p      -iat 

-2:tia    /  e-^^^{i  +  |(t-v)-  (t^-   tv+v2)}dt-  2n±c.v\^f    I     f-^^  dt 


(42)      -   2rtia    /      e"^^M|+|(t-v)-(t2-tv+v2)    +-^    (^)Hdt 


.2«(-_^_v    )    4---    (-  +  V)    --   -   2:tiav    (--)    j     

o 


t+v 


dt 


-   2rtiae 


3  "^1  / 

1."^-^^    r  t^(i-it)^/^-^v 

,/         1+  V-  it 


on  evaluating  the  first   integral  and  deforming  the   contour  of  the  third  integral 
into  the  line  joining  1  to  1   -   i  <».      Hence 


1 


^   t    ^  t+v 


e  dt   =  ^   (a^v'^+  iav-  2){^^)^    -   2rtiav-'(^^^) 


X 


-iat 
e 

t+v 


a 


dt    -   2jtae 


V 


00  1 


5 

1^^"^°-  r  t^d-itf   -at,, 

/      ; — ^ TT—  e       dt 

,/       1+  V-  it 


so  that   P-^(v)   satisfies 


P   (v)      =     2ae 


1    •    •       .  "'I  ^ 

T5Ti-ia-iav  1       p    ,2/-,     .,  n2  , 

-X_^2       /     t^(l-it)        -at  _    /_v_)-  ^_ 

1+ V-  it      ^  ^        '^l+v^  Jt 


(— -)' 
^1+v^ 


V    v'^   e 


-lav 


1  -D   f4.\    -io-'t 

i     Pn(t)e 


X  MV) 


l-t^t  1 


t+v 


dt    . 
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On  putting  P-,(v)   =  (- — )^  e    R(v)  we  obtain  an  equation  analogous  to 


(31)  of  Jones  [17],  namely 


1  .  .    «>   1       o 


n/  ^    o   ^      /  t^n-it)"   -at,^    1   r  ,l-t.^R(t)e 

R(v)  =  2ae       /  r-f -^  e   dt  -  -   /   (777)     ,  '   dt 

.  /   1+  V  -  it  Jt  ,  /    1+t     t  +  V 


o 


Therefore  we  form  the  iteration  scheme 


0°  a  /^  \  -2ia(t-v ) 


R  (v)   =   (^^=-^)''   /  -^^-— dt  , 

n^  '      ^  J  t  -  V         ' 

o 


iay  p  X 

(^3)      Sjy)   =  ^^tr"  j   (^)2  Sn(x)  M(x,y)dx      (n>l), 

o 


r  ,  ,  i  s   rt)e-2i^* 

>^)      sjv)  =    i^f     ^^^^^ dt  (n>2). 


(U5)      s^(v)  =  2ae       /  ^^^  .^^  e    dt 


On  accoimt  of  (^0)  we  have 

n  w^(l-w)2  ^^  \a  ^^ 

'  /  N  io-w      .        00     T    n    fj.  \    -2iat 
S  (w)e         law      r  t-1  i  S  (t)e 

n=l   L  1 
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Now 


1 


2_      /    (lzV)2  g^v^-lavrg  ^^ 


1 


V    -   w 


a 


Q        ;     1  -A7  ^   I    2        ,  N     g  w  -iaw-2V 

2      /    ^V^^     <;a  v-a(i-aw)  + 77^r^>  ^^ 

a      ^ 


2rt 


i  a^    -  Ja(i-aw)  -  (aV-iaw-2) 


and 


1      o-™,,-* 


,1-Vst  2iav- 


•3,       o      -iat 


^    V  V  -  w 


t    -   V 


dt   dv   =  -2ia        e 


■iat    /    /1-v  \2 


2,,      ,    2 


X      /v    +   t    +  w    +  (^_w)(v-t) 


t   +tv+w    )-\rt(t+v)y     ^^  ^^ 


=  -2ia        e 


■iat 


O  2 


1  1 


-n+«M(l-w)^5(t-w)+«t2(t-l)^   ^^^^^ 


w 


wt(t+w)^«    (^^)^&(t-w)+rt(^)        t 


w 


dt 


from   (e6)  and   (E3)   of  Jones    [iTl-      The  last   integral  can  be  written 
-2«ia    re"^^N^+|(t+w)  -ft2+tw  +  w2)>dt    -2ia«V   e""™(l-w)2 


-2«iaj     e''^'^<(|+|(t+w)-(t 


^.tw.w2)-T!^.    (¥)'W^- 


t-w       t 


2h  - 


A  comparison  with  (1+1)  and  {h2)   reveals  that  we  have  the  same  integrals  with 

-w  in  place  of  v  and  since  no  principal  values  are  involved  we  can  deduce  that 

5 

2  ^  axe^^lllillZi  r  (l-it)Ve-"^  ,, 

p(w  j  '^^  2iaw  +  i ; /  -^-i; ^— n dt 

i/i   ^4     J  1  -  w  -  it 
rtw^fl-wj'i      "^ 
o 


(U6) 


„  ,  ,  -iaw  .        00       a    /^\    -2iat   ^ 

S  (w)e  law  p     i  S  (t)e 

+  ;    <  -S +  -^ r  /  {^y  —^ dt 

w  ^/,   ^?  J    "t      t  -  w 

n=l   I  :rw^(l-w)2  ^^ 


With  p  known  p  can  be  found  from  (38)  and  the  whole  distribution  is 

known  as  soon  as  the  constants  d  .d-,  and  A^  have  been  determined.   We  turn  to 

o'  1      1 

this  problem  in  the  next  section.   Before  doing  that,  however,  we  quote  the 
following 

.   /    -,  ^    1   .       00  11 

Pae^"^^^"^^  '^"^^    r  fl-it)2t2    _at 
f(w)~2ia  +  ^^^— T i /   V"     •!   e  °-^  dt 

2/n    V2        J    1  -  W  -  it 

nw'^(l-w)'^      ^ 


{hi: 


,,  /  %  -iaw      .        00       TT  /i  ^  -2iat 
H  (w)e  law      p  .  -,  i  H  (t)e 

^ -1--^ r  [{^f^ dt 


where  H  and  h  are  related  by  equations  of  the  type  (^3)-  (^5)  and 

1   .     00  11 

h^(v)   =  2ae        J   l+^.'i^    e    dt  . 

o 
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5-   The  equations  for  d-  , d  and  A.^  . 

It  is  clear  from  (3^),    (i|6)  and  (kl)   that  d(w),  p(w)  and  f(w)  all 
become  singular  like  (1-w)  ^  as  w  -» 1.   Suppose  that,  near  w  =  1, 

f(w)  =  {l-v)~^{->^+X^(l-v)+  X^{l-wf]  +  Oiil-vf]    , 

d(w)   =   (l-w)"2(^^+  H2(l-w)+  ^^(l-w)^}  +0[(l-w)^}, 

p(w)   =  (l-v)'^[v^+  v^{l-v)  +  v^il-wf]    +0{(l-v)^}  . 

It  is  not  perhaps  immediately  obvious  from  (3^)>  (^)  and  (kl)   that  there  are 
no  terms  involving  the  powers  (l-w)  and  (l-w)  but  this  is  proved  in  Appendix  C. 

We  have  seen  at  the  beginning  of  Section  k   that  p  (w)  must  vanish  like 

3/2 
(l-w)  '   as  w  ^  1.   Therefore  (38)  implies  that 


[hS] 


q  i 
d^X^  -  A^n^/(8m^)2  +  d^v^  =  0  , 


'k9)   d^X^  -  A^^^/^S^^)^  +  ^1^2  "   ° 


Also,  it  was  pointed  out  in  Section  1  that,  because  of  the  edge 
conditions,  a-^^(r-j^)  +  h-j^(r-|_)  =  0  at  r-j^  =  1.   On  account  of  (29)  and  (36)  this 
implies  that 

-^   {A  d(w)  -  f(w)]  +w2(-|-)  p  (w)   =  0 
1+m^    1  wdw    1 

at  w  =  1,   This  requires  that 
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(5o; 


0  . 


Eq.uations  (1(6)  -  (50)  constitute  three  linear  equations  for  the  constants 

d  ,d  and  A  and  hence  determine  them  (so  long  as  the  equations  are  independent! 


It  is  shown  in  Appendix  C  that 


\ 


2a  h 


\     -^ — j^  •'      e   dt  + 
^    t^ 


n=l 


00     ,  >  -2iat 

la   p  H  t)e 

V  \  -^^ r-dt, 

"    ^.  t2(t-l)2 


^2  =  ^2-^''K^~   ^ 


K-,(t)     ,       1 
r~  e   dt-  23r  a 


t' 


^e 


la 


dt    Pn(^)^ 


-2iatl' 


n=l 


:t-l)' 


t2 


A   1.    1  2>     .1   .  ^,     2a 
^3  "  ^8  ?°  ?  )X^+  (-  -ia)A,2  +—  e 


t-Li,!   r  Millie  .-t 


-jT^a  - 


e   dt 


2 


la 


n-1      1 


dt 


(t-l)2    t2 


H  (t)    ^.  ,  , 

n     -2iatl 

e 


2e 


la 


^1  =   « 


CO        ,    V 

ia   p  E  (t)e 
e     /   n' 


-2iat 


n       J  !_  L 

n=l      1   t2(t-l)^ 


dt  , 


=  (-  -ia)M.^ 


2e 


la 


dt 


E  (t)e 
'  n 


-2iatl 


n=l 


:t-i)^ 


t2 


> 


,11.    12,     ,1      .     .        , 


he 

3ir 


la 


dt 


'^n^*^   -2iat 
~r~   e 


n=l 


;t-i)2     t' 


]■ 
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1    00       ■^              .    °°  ^  / ,  N  -2iat 

2a  U   /  (1-it  -o-t,,  ^  \  e /  n\ ,, 

V-,  =  —  e    /  -^ — f— ^  e   dt  +  >  -—   /  T  T^  at, 

1     «        J      4-2  Z_l   «  J  +2/4._-,')2 

o    *  n=l  1  ^  ^*  ' 


2e"°^   r   dt 


n 


n=l   '   i   (^-l)"!   ** 


where  K  (t )  =  [(l-it)  ^-   l}/t.   Since  {k8)   is  satisfied  the  terms  involving 
X,  ,    \i     and  V  in  A.p,i-i  and  v„  respectively  can  be  ignored  for  (^9)'     Because 

{k8)   and  {k9)   are  satisfied,  the  tenns  involving  A^  ,  Xp , n-,  , |.ip ,  v.,  and  v„  in 
X^jH^jV^  can  be  neglected  as  far  as  (50)  is  concerned. 

The  above  formulas  permit  the  calculation  of  A_  ,  Xp, .  .  . ,  v^  to  any- 
desired  degree  of  accuracy.   The  first  few  terms  of  the  asymptotic  development 
are  derived  in  Appendix  D.   Thus 
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(51)  \ 


2a'=^e 


1    . 


32a        128a- 


2ia  -  7rti 


+   0(a'3)    , 


la 


(52)   ^,    = 


1  It 


2e        /-,      i_ 1      N   ,  e 


Titi-ia 
4 


128a 


'J 


■n     a' 


-2ia  -  Trti 
512a        l6n  a  ' 


1    . 
1   Trti 


-2ia  -  rni 
-2ia             I ,  .                       k  , 

K.)    ,    e /-,  _  itli  +  e ^        oCa"-^) 


(53)v,=S^(l-a.-a-, 


(5^)    Xg   =   (2  -ia)X^ 


3       1    . 
2     -Titi 


^2  ^     32a^     128a-^ 


-2ia 


■1+^+^       1 


-2ia  -  73ti 
4 


16^  ~7ZU/^ 


l6jt2 


•)+0(a-^). 


a" 


-2ia  -Trti 


512a  Ibrt  a  ' 


,^^^  /I      .     N  le ,  4ia 

,55)  ^2  =  (i-^'^V  U^^^'JT^ 


i        1 
2     -TJti 


.    -2ia 


(56)    V. 


,1      .     ^  ija       "1|      ..,   ,  5i   I     15         Hi    N   _^  ie /,      iili.e 

(2-xa)v^--^e  d  ^  t  "^r^  "  "irS  ^    "        n  ^^'iG^'^.A 


-2ia  -  Tjti 


32a        128a- 


I6rt2 


a' 


372" 


^/    -2n 
+  0(a      ), 


(57)    X,    =    (i+|ia  +  |a2)A^+  (|-ia)X2 


3 


ijo-e 


1  1    ■ 

2  I''^ 

3^2  '^     32a^     128a^ 

-2ia  -^Jti 

-2ia       ,^      5i        e n    ,    ^/    -1\ 

. (1-T^+ i  o  /o  '    )  +  o(^      )^ 


3n 


I6a     ,.k3/2 


l6rt2a- 


(58)  M 


3 


(59)  V 


,11.        12,  ,1      .     ,  I6a3/^     t'^^-^" 

(q+21^+^    )u^   +    (--  la)^2   -  :^JJ^  e 

,11.        12,       ^,1      .     ,  8a^      h'^^i^^l^       45    ^    15i 

(^+-ia  +  ^    )v     +(--ia)v^--T  e        (1+-^;-— ^+-^ 


-21a  -  Trti 
^^     512a^     I6rt2a^/^ 


3         ^8      2  2      '1      '2 


3«' 


^     32a2  '^gSa^' 


hi.     -2ia       /^ 
-^  e  (1 

3rt  ^ 


-2ia-TJti 
^  +  ^-^ T—)    +   O(a-l) 
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Substitute  (51)  -  (56)  in  (¥3)  and  (^9)-   A  simpler  equation  is 
obtained  by  forming  the  combination  (hS)   +   {l/kci){k9) .      Then,  if  (k9)   and 
(kS)   +   (±/kx){k9)   be  solved  for  d  and  d  we  find 


^       12&x2 


1 


-rnl-3±a] 
512a     Ibjt'^a  ' 


J 


I   1   1  •^• 


(61)   d^(i- ---^)  =— -^<;— e         (i  +  &+7;:t:2^ 


6iKx     U(8m^)2   jr2  ^  128a 


le ,   5i    127  N   e 

jt        ^        l6a      _^  2^    "   ,.  1  3/2( 
512a  Ibn^a  '     ' 


Next   subtract   i|ia/3  times    {k9)   from  the  term  in    (      }   in   (50)  and  use 
(5^) -(59).      Then  substitute  for  d     and  d     from   (6o)   and   (6l).      The  result   is 

jt     1  I         Sx     2^i^i  l6aY  ^2  ^'  i^no. 


Hence 


1    .    .       r  TTti-2ia  i    .    ,  -  fni-6ia> 


-I    i  i  TJtl-ia    I  4^  ,  -g-nl-i^-la        ^4  -„ 


l^"^""^"  ,,.  irri-l^ia  V 

P^2„2  l6a  k-ncL  .    o   o 

2«  ^  Bit  a'^ 
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Equations  (60)  and  (6l)  can  now  be  used  to  determine  d  and  d  but 
these  constants  will  not  be  required  for  the  subsequent  analysis,  although 
they  are  necessary  to  specify  the  distribution  on  the  disc  completely  via  (38 ^ 


6.  The  scattering  coefficient. 

It  is  known  (see,  for  example,  Jones  [23]  )  that,  if  an  arbitrary 

-ikRnQ-n 
obstacle  irradiated  by  the  plane  wave  with  E  =  e  e        produces  a 

scattered  electric  intensity  E  ~^  (n)e    /r  in  the  direction  of  the  unit 

s     s 

vector  n,  the  scattering  coefficient  is  -  77;-^  A  (n  )'e  where  S  is  the 
-■ '  *  kS  ^  -s^-o  -o        o 

o 

area  of  the  incident  wavefront  intercepted  by  the  obstacle.   Translated  to 

our  coordinates  in  which  the  radius  of  the  disc  is  unity  this  states  that,  if 

-iaRriQ-n  -iaR/ 

the  plane  wave  E  =  e  e         produces  the  scattered  field  E  ~  A  (n)e    /R 

the  scattering  coefficient  of  the  circular  disc  is  -  — ^  A  (n  )-e  . 
^  a  '^  — s  o  "o 


I.e. 


Now  it  follows  from  the  representation  (l)  that 


P  -iaR  n      ian-R 


Ag(n)  =  -an  a  <n  a  /  J(R)e       dS) 


But 


n  =   (sinG  ,O,cos0  ),  e  =  (cos  cp  cos  0  ,  sin  cp  ,  -  cos  cp  sin  0  ) 
— o        o       0—0        o     o       o        o      o 


and  so 


1  23T 

A(n)'e     =a      /      /{(Jcoscp-Jsincp)coscpcos0+(Jsincp+Jcoscp)  sin  cp    } 
— s  — o     — o  J     J  ^  9  o  or  *^        cp        ^  o 

°     °  -iar^  sin  0   cos  cp 

1  o 

X     e  r  dr  dcp   . 

11^ 
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Substitute  for  J  and  J  from  (3)  and  use  the  formulae 
r      CD 


27t 


imcp+io-r   sinG  cos  cp 


1 


dcp  =  2jte  J   (ar,  sinG  ), 

m       1  o 


Jv-l(")+  Jv+l^")   =  2vJ^(2)A.    J,.l(z)   -  Vl^^)   =  2j;(z), 


Then 

00  1 

„    /       ^  2  V^     i(n+l)m 

A   (n    )-e     =  Jta       >      e^ 
-s  -o     -o 


n=-oo 


nJ   (ar^sinG 
n       1  o 


;a   -b    ,  .     ^ 

n     n         ar^smG 
1  o 


[a  +h    ) J' (ar^  sin  G    )  >  cos  G  cos 
'   n     n     n^      1  of  o 


(63; 


-> 


nJ   (ar   sin  9    ) 

i  <  (a   +b    )  — ■      „     °     -    (a    -  b    )J'  (ar,  sin  Q)}  sin  T' 

^^   n     n         ar   sm  G  ^   n       n     n       1  of  o 


r,  dr. 


from  which  the   scattering   coefficient   can  be  determined. 

In  the  particular   case   of  the  normally  incident  plane  wave  G     =  (p     =0 
and   (63)  reduces  to 


A   (n    )'e     =  ivx        /    (b^  +a  .,  )r  dr^    . 
-s'-o     -o  J        1        -111 

o 

Since  a   =  b  the  scattering  coefficient  o   is  given  by 

1 
o     -     -8jta  ^  /  b  (w)wdw 
o 

or,  from  (29), 


(64) 


=  1^  /   (f(w)  -  A^d(w): 


wdw 
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7.  The  scattering  coefficient  at  high  frequencies. 

1 

The  quantity  (2/a)<^  /  wf(w)dw  is  the  scattering  coefficient  of  the 


sound-soft  disc  in  a  normally  incident  plane  wave.   Its  value  at  high  frequencies 

has  already  been  calculated  (Jones  [l?]  )  and  it  may  be  denoted  by  a    .   Thus 

1 
(65)  a     =     a  -  -^  A        /  wd(w)dw  . 

o 

For  a»  1,  d(w)  is  given  by  (3^).   The  first  term  in  the  expansion  of  d(w) 

gives  a  contribution  to  the  integral  in  a  of  a  where 

1    1  . 

I      /   2w  e       -, 
a       =   / dw 

°    J      jr(l-w)2 


(66) 


ia  -  jn±  ^rti      _  5 


The  general  term  of  the  series  for  d(w)  gives  a  contribution  to  cr  of  -  -q)  A-,  a 

where 

If  1  .      <»       T7  /+  \  -2iat 

p         .      2  lo-w    p  +  n  1  E  (tje 

a'   =  He  (w)e-^"^  +  ^^^— r  f  i^?  —, <it>  dw  . 


The  asymptotic  behavior  of  cr  may  be  determined  by  the  same  procedure  as  was 
used  for  the  corresponding  quantity  in  the  scalar  problem.   All  that  is 
necessary  to  put  e  for  h  .   Consequently 
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-la 


2  . 


,(1)+T&e   (D+^e'd)-!:^ 


n  2rta      |      n^"^^'^l6a  ^n''"^'''^^  '^n'' 


e^(0] 


-V^ 


(67, 


+    ■  i  3/2 


n^    ^     a   (16     n 


(0) 


2  <(0! 


-jTi+ia 


^^ 


8rt' 


3/2 


e      0 
n 


tf„f°'*?<(°>' 


+  o(-t; 


a 


Because   of    (DI6), 


■la 


a     =  ^= 

1  rta 


+  IL__2li.        e_ 


-l^^l 


1       .        . 

Trti+ia 


i6a  -  &  ^  -^  XsT^ "  7X37^ 


+   0(a 


5 

2, 


rt^a  '  ij-rt^a- 


and,    from   (Dll), 


(69; 


le 


-la 


Srra 


1   - 


^       -i 

2^)    +   0(a  ^; 


(70)        a'      =     0(a-^; 
n  ^        ' 


(n>  1) 


Combining   (66),    (68)   -    (70)  ve  have 


wd(w)dw 


2e 


la  -  -Tjti 


1    1 


'^-k' 


-la 


jta 


l-?fc)    +  0(a2) 


16a 


Hence,  from  (63)  and  (65), 


-  3^  - 

1  3 

7jri-2ia  i    .    ,  .  ,  rtl-6ia 


"  =  -  ^.  -!S<|l-^^^TT-(l-*)-S*=----  -   " 


soft  "a'^V   ^ii   ^   ~  l6a^  "8a     l+rox      „  3/2  3/21 


TJti-2ia 


Since 


cos(2a-  TJt)|       _-, 


it  follows  that 

0-  =  2  <  1-  I    o  /o   sin(2a-  TTr)  +  — 5+ 5  sin(lw  -tjt) 

1    jr2a^/^         ^    Ija   2jta  ^ 

(71) 

1           3      27  1   I      - 

"  — n/p  p:/p  sin(6a-^n)-  1  s/p  cos(2a-  -rjt)  >  +  0(a 

W/  a^/^         ^    I6jt2a^/^         ^  J 

This  scattering  coefficient  agrees  with  that  found  by  Seshradi  and  Wu  [15] 
except  for  the  last  term.   The  numerical  factor  of  the  last  term  is  foimd 
by  them  to  be  T/k-   instead  of  2j/l6;   the  difference  is  probably  due  to  an 
arithmetical  slip  on  one  of  our  parts  but  this  has  not,  so  far,  been  located. 

Since  all  the  high-frequency  integral  equations  solved  here  are  of  the 
type  for  which  uniqueness  and  an  asymptotic  development  were  proved,  the  same 
properties  are  valid  here.   Consequently,  it  would  be  possible  to  make  an 
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estiraate  of  the  magnitude  of  the  order  term  in  (71)  thoiigh  it  would  be  rather 
complicated  because  of  the  intervention  of  d  ,  d  and  A  .   Experience  indicates 
that  a  formula  such  as  (71)  may  be  used  for  numerical  values  when  a,  >  3  and  will 
give  tolerably  accurate  results  even  for  an  a  as  low  as  2. 


APPENDIX  A 


In  this  Appendix  is  derived  the  formula  for 


2.2^         _  _      ..  ^2 


2rt     -ia.(r    +  ^r    -   2rr   cos  cp 

r>  11 


t   (r,r    )   =      / 2 2 T~  <^°^  ^^^ 


r    +  r,  -  2rr,  cos  cp)2 
u  11 

We  deal  firstly  with  the   imaginary  past.      Now 

2Tr      .     !    ,    2        2      ^  ^^1  2rt    _.,  r    /    2  ,      2  ,2't 

P    sinla(r    +r     -  2rr   cos   cpj    J  -,        i     p    Jila-(r    +  r    -  2rr   cos   rp)    J 

/ t; cos  nrpdcp  =   (^ir)^    /     ~ ~ Y —  cos  ncpdcp 

^  (r    +r,  -  2rr,  cos  cp)^  ^  (r~+  r,  -  2rr^cos  cp)2 

oil  o  ^  1  1  ^' 

2jt     i^ 
=  a    /      /   sin  0  J   (a(r    +  r,  -  2rr  cos  fp)"sin  0)   cos  ncp    dG  d9 
o     o 

after  using  Sonine's   first   finite  integral   for  Bessel  functions    (Watson    [2U]). 

The  addition  theorem  for   Bessel  functions  gives 

i_ 

2Tt       2rt  0° 

-^if     =  o.  /   sin  Qy      e  J   (ar   sin  0)   J   (ar     sin  0)   cos  mcp  cos   ncpdG  dcp 

o     o  m=0 

where    e     =1,    e     =2(m/o).      Hence 
o  m  \      /       / 
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1 

23t 


•^  if     =  2jta    /      J   (ar   sin   0)    J   (ar^  sin  0)    sin  0 
n  .  /        n  n       1 


dG 


Since 


^zdz'^  V  ^       -^  v+m        '         ^zdz'^  v  v+m 

z  2  z  ^ 


and  therefore 


1 
2n 


\n  n  __^_.'^  In  [   ^   ^.^  Q^  j  ^cLr,sin0)  sin^'^^G 


■^11/   =  2it(-)  r"(-~— )  ct^  ''   J  (ar  sin  ©)  J  (ar^sin©)  sin  "&i0. 
^  ^n  l^r  dr        J   "^  o   1 

o 


It  has  been  shown  elsewhere  (Jones  [22])  that 


1 
2rt 


J     i(z  sin  0)J   i(Z  sin  0)sin^"^0  cos^^G  d© 


o 


i«  T    f^„2_  „2^^  2  .2 


lip  J  ,^Uz  -  Z  cos  ©)  ,      ,  -,      p-x 

=  z-^Z--/  J^_,(Z  Sin  e,   T  ,      .(.,,)   =in"-^^^«  eos=^«e 

^  (z  -  Z  cos  0/ 

o 

when^(X)  >  4,  %i^)     >  %(^)   -1,  '^(^)  >  -2.   Use  this  formiaa  with  X=0,   U  =  n  -  i, 
V  =2,    Z  =  ar,  z  =  ar-,  and  then 


-■  2   2   2^.|. 

-±.   P  ^.^L--vr^  -  r  cos  0)"  J 


1   p  Jl{a(r  -  r  cos  9)^j      x 

-^.^  =  2.(-)V;(^-^ra--\ar)^  /  J^_^,(ar  sin  ©)   \  g   ^   2,,l/4  -""  '  ^^^® 


■1°-^1  ^      "  a,2(r^  -  r  cos  0! 

o  J- 


fAl 


^n-l^n  V^dr^    J   n-2  (r^-x^)^/^ 

o  -L 


on  putting  r  cos  ©  =  x. 
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Turning  now  to  the  real  part  we   employ  the  result    (Watson    [2k] ), 

P    J   (bt)t 
cos  by  /        o J, 

and  then 

2jt  00 

r  r  2       2  -              t 

%ii     =  /  J    fr    +  r!:-  2rr   cos  9)^t)   — ;- — r-f  cos   ncp  dt   dcp 

n      J  J  o              ±            1  (f'-a    )2 

o  a  ^ 


2rt       00 

t  / 


7      e   J   (tr)j   (tr^  )cos  mgi  cos   nq)  dcp  dt 
/    I    mm  ml  -^      -r 


(X      ^  '        o       m=0 


2n    J     -^^J,(tr)J^(tr^)dt 

^        (t    -a    )^ 
a 


l^r^dr    '      ,/       f.2     2  ^  2     n^  o^      1 

11        ^       (t   -a    )■" 
a      ^  ' 


2"^-)  ^i^r^)   J  ~rjjr-  ^o^^^  ^  )  ^i^  ^* 

11  (  t     tCO     j 

o        ^   ■ 


According  to  Sonine's  second  finite  integral  (Watson  [2U] ) 

2  2  -         - 

J  {(t  -kx  )^r)    ^  r^"  J  i(tr  sin  0)  J  1  (ar  cos  G)    1       1 

n   n  /    -2'  n-2             p            n+^ 

o     TTZ —  =  r   /   — '^ 1 '^ 1 sin^G  cos   ^GdG 

(t^^a^p       J      (tr)-2  (arf"2 

and  hence 
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%\  = 


1 

Ini-)   r^  ( — -; — j   ra  /    J      i  (ar   cos   Gjsin  0  cos        0X 

o 


n-2 


t'J  j_(tr    sin 


»)J    [(t^-Kx^)2r    )dt   d0. 


Now 


r 


—  ?     ?  — 

t^Jl(bt)j    {a(t    +a    )2}dt 


a  <  b. 


a  >b 


and  therefore 


2      2.2, 


^  n  1  r  dr 


J  i{a(r   -r   sin  0)"^} 


n  n,      d      \        2   1-n        -,        ,  _s        n+2'       -^ 

"""       "        "   '  J     1  (ar   cos   0JCOS        0  ^ — ^- ^ — -,  rr 


d0 


(r   -r   sin  0, 


the  upper  limit   of  integration  being  ~n  if  r,    >  r  and  sin      (r   /r)   if  r     <  r. 
The   substitution  r   sin  0  =  x  gives 


(A2: 


H^-^-i-f-^n^ 


n     "1 ^ d_     \n 

r^  di 


min(r,r    ) 


1    J  ^(a(r   -X    )^ } 


a        r 


Y'l 


J     i.{a(r"-x    r] 
n-2 


2      2x^.-2      '1  "    '    '    iJ2     2^t(n-2 

2  , 

r^    -   X 


1 


r-_x-)2(n-2)^^ 


Combining    (Al )  and   (A2 )  we  obtain 


*n 
(A3 


n 
r 

=  2jt(-)     :; ( — - — ) 

n-1  n  V  dr^ 
a        r  11 


n 


mm  r,r-,  j  ^     „   i  /      i  x 

r         T        f    r    2     2.i.,      .    ,    2      2i,    (r^-^2)^('^-2) 

o  ^   1        ^ 


,2      2,i,    fr^-x^)^(n4) 


i    /   Jj^_i{a(r   -X    )2}Ji{a(r^-x^)^}    '^    g^^  ^  j^y|^       dx 


(  i^3_-x    ) 
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Since  ^1^  is  symmetric  in  r  and  r  we  can  interchange  r  and  r  on  the  right-hand 


side  of  (A3)  and  it  will  still  hold,  i.e.. 


n     ^     n-1  n  ^rdr 


(A4) 


a   r 


1 


.2  ^2.i(n4) 
J      J,..(-(ri-x-)^]J_l(a(r2-x2)*]  ^J./,       dx 
u  o  (r  -X  )  ' 


min(r,r  ) 

r  J  t--2„2.i 


""l  /  2  2,i(n-i) 

-  i    J  i{a(rf-x^)^}Ji{a(r^-x^)2]   ^^   ^  w,    dx 
,  /    n-2    1       2 


/  2  2.1/4 
(r  -X  )  ' 


APPENDIX  B 

In  this  Appendix  we  indicate  another  approach  to  (15)  from  that  given  in 
the  main  text.   Since,  for  n  >  1, 


,  n  n  /  2  2vn-l    ,    .  ^  ^     n    ,    .,   2     2vn-2    ,    s 

{   (n-1);  2^-1   u^-1        ^^^    {   (n-2):2'^-2    ^n-1 


(-^) 
Vdr^ 


4-  /  i-f'^^du 
rdr  ,/  n-1 

o 


(-)   g(r)/r 


(15)  is  equivalent  to 


^0 


"■  J_i_(a(r2-x^)2}   ■  ^^  Ji{a(r2-x2)2} 

'  2   2,lA   ^(^^^^  -  ij   ',  2  2,lA  ^(^^'^^ 
Q    (r  -X  )  '  ^     (r  -X  )  ^ 


(Bl) 


^^^    o^       n-1  p  ,  2  2.n-l 

„  2m      g /   u  -r  )     /  s  ., 

C  r   -  /  -^ i^ g(u  du 

_          n-l);2  jt  "^      u 

m=0  o 


where  C  ,C-,,...,C  .,  are  arbitrary  constants.   If  n  =  0,  the  right  hand  side 
o  1      n-1  ''  ' 

is  replaced  by  ag(r)/2jt. 

Replace  the  Bessel  functions  by  their  trigonometric  equivalents  and 
then  (Bl)  becomes 


l__.._r_  /  2  2.^- 


cosfgfr^-x^)^]  p sin{a(r^-x^ )A 

2     2- "^M^^   -   i  /  2 — 2~~^ h(x)dx 

'-^    (r  -X  )^  ^     (r  -X 

o  o 


n=l  1        r   ^  ^   ^ 

„       n-2       p  ,  2  2  sn-1 

(n-l):2^^V-'  ^     u"-l 
m=0  o 

where  C  =  (rta/2)^C  .   Make  the  change  of  variables  r  =l-z,   x  =l-t,  with 
mm  ' 

the  result 


1  i  1  i 

^  cos(a(t-z)^}   (^^^^  _  .  r  sin[a(t-z)^)   ^^^^^ 

(t-z)2     1         ^^    (t-z)^      1 


o 
(B2 


n-1  ^  i        1       ^ 

C'(l-zf-   ""'  A^   r  -^^^^^g((l-u)*}du 


kl 


where  h  (t )  =  h{  (l-t  )'2"]/2(l-t  )2. 


Equation  (B2)  is  of  the  same  type  as  (2^1-)  and  so 


(B3: 


^^1^^  )  i        w+v 


-ia(w-v ) 


w-  V 


dw  =  H^°\v 


where 


h'°'(v: 


1  ^ 

2  dv 


n  2    ' 

/   cos(a(v  -z 


2   .t- 


,  2 
(v  -z 


1-1 


C'(l-z) 
m 


m 


a 


1 
n-2 


m=0 


(n-l):2 


n-t^  i 


.n-1      j^ 
"""^   -g{(l-u)2)du 


:i-u) 


1 
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dz. 


The  behavior  of  h  ,  and  hence  of  h,  for  large  a,  may  now  he  determined. 

The  constants  C  have  to  he  chosen  so  that  h(x)  ->  0  like  (1-x)  as  x  -* 1 

because  of  {l6) ,   when  it  is  remembered  that  k(r-,  )  is  nearly  a  constant  multiple 

of  (l-r^  )  ^  as  r^  ->  1.   Thus  the  C  are  chosen  so  that  h,  (t)  -^0  like  t  as  t  -*  0. 
11  m  1 

Over  h  is  known  k  is  found  from  (l6)  as  follows.   Since 


-^  {z"j  (z))  =z^-lj   (z' 
zdz     v^  v-1^    ■ 


d 
xdx 


f   2  2.i(n4) 


J  i{a(r^x^)2} 


-a(r^-x  )''     J   {a(r^-x  )'^) 
"2 


Henc« 


1(2 


(S;'   -(>=) 


/    ^n  n    r  ^'    .    .    2      2,5-,    /    2     ; 
=      (-)  CO      /    J_i{a(r^-x    )'']    (r^-x 


2..2,i,    ,^2_^2,-lA^^^^)^l-n^^^ 


i_       p    cos{a(r^-x    )^}     k(] 


(_)n^n-t   ^2^t      /■    ^^"^    ^^1  "    '    '     irr 


p     p   1 


n-1 


ir. 


from  which  can  be   deduced  that 


(B4) 


n-1 


f-r^  .2,i  d 


e)' 


n-2   rt   dr^ 
a   '^         1 


1      v,f  /  2  2,i 
n  X  coshia(x  -r  j  , 


(^)  h(x)dx 


An  alternative  way  of  writing  this  can  be  obtained  from  the  fact  that 


n  X  cosh{a(x  -r^)^}   ^j  ^  1 
r.dr.  ,/     ._2  __^A  ^^^   ^^^^^^ 


11 


(x  -r^)' 


1   _d_ 
r.,   dr., 


-il 


a  "^"^  ^  {a(x2-r2)2}(-|-f-lh(x; 


1   ^_  Tx. ....,_  r_. /__2  __2,i 


■^ —  /  -sinh{a(x  -r^  ) 
r  dr  J  a  1 ' 


X  cosh[a(x  -r  )^ } 


1^  '    ,  d  s^i    .^ 
2   2.i (^^   h^^)^ 


(x  -r^)' 


since  h(x)  has  a  zero  of  order  n  at  x  =  1.   Hence 


(B5)  k(r  )  =  V^  ^?)'^" 


1  2   P   ■ 

/  nH+I    1          ,T  p  X  cosh  {a  (x  -  r^^i^' 
(- )     ^2^2  n,   d   ^n+l  /  ^_   _!: 


n-2    ^jr   1  r^dr^  ' 
o.  11 


.2  2>t 
(x  -r,  ) 


h(x  dx 
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APPENDIX  C 
This  Appendix  is  concerned  with  evaluating  the  behavior  of  f(w),  p('w) 

and  d(w)  near  w  =  1.   It  is  apparent  from  (kj)   that  there  are  two  different 

types  of  term  to  be  considered  for  f,   namely  a  typical  member  of  the  series 

and  the  member  outside  the  series.   There  is  a  distinct  similarity  between 

{k6)   and  (hi)   and  we  can  deal  with  the  first  members  of  both  at  the  same  time 

by  evaluating 

ia(w-l)-77ri   <»       i   i    . 
2i,,i(m-l)  ,  2a^_ _   r  (l-.t)^te    ^^ 

o 

and  then  giving  m  the  values  1  and  5  respectively  for  f  and  p. 
Now 

f  (l-it)Hi-"'^  ,,      .  f    (l-it)^  -at^^    .  .^  ,    f  e-^^     . , 

/  ^— '■ — dt   =  1  /   ^ 1  —  e   dt  -  i(l-w)  /   "T         dt 


1-  w  -  it 


t2  ^   t2(l-w-it) 

o  o 


p  K  (t)    ^  _  p   1  .  -at 

+  (l-w;  /  1—  e   dt  +  (1-w)   /  -^ dt 

J  t^  ^   t2(l-w-it) 


00       /    \       1 

„   p   K  ( t  )  +  2^1      4- 

(1-w)   /   — Y s    clt 

^    t2(l-w-it) 


1 
2m 


where  Kjj^(t)  =  {(l-it)"^  -  l)/t.   Also 


kk 


-at 


t^fl-w-it' 


dt    = 


ijt  e 


2_ia.('l-w)     r    e 


-i(l-w)x 


dx 


a 


1    .    /-,       N       -  Tni    1  1 

.      2     ia(l-wl      J         k  2/n  \-2 

lit  e  <  e  JT    (1-w) 


af 1-w) 


■IX 


'1-w) 2x2 


dx 


=  Tte 


Tni+ia(l-w) 


1-w: 


-i-^ 


3 


1 


2o}-^<x{l-v), 


+  0{(l-w)^} 


Hence 


fl-it)H^e" 
1-w-it 


,  ,  N2m 


dt 


-     r    (1-it)^'"      -at,,         ,,       , 
i    /      "^ 1     —  e        dt   +   (1-w) 


t2 


K   (t) 

m  -at,, 

— T —  e       dt 

t^ 


*-^  /        \  1 

2    r   \'.^)+2tni      _^^ 
-   i(l-w)    j     —37^ e        dt 

o 


+  <rmi(l-w)    -  i(l-w)  >  Jte 


) 

(l-w)"2.  ^T-    {2a2-^a2(l-w)} 
^2  J 


+  0((l-w)5/2) 


Therefore 


h^ 


2  law 


2(m-l)    ,      2ae 

+  1 


ia(w-l)-  Tj^JTi 


rtw2(l-w)2 


.  1-it)      t   e 
1-w-it 


dt 


2ae 


V 


TTl 


Tt(l-W) 


1  ^11'^  2( 

,1+  (--ia)(l-w)  +  (^  --ia  -ra,^)(l-w) 


'CI 


' 1  —  e       dt   +   (1-w) 


K    (t) 

T—  e       dt 

t" 


^    P  K  (t   +:3ii        , 

i(l-w)      /      T77^ e        dt, 


t 


JW 


i   1. 

,2^2 


2ia 


-   jt^^^^l-w)  <2  -  (^+m  -l)(l-w) 


Turn  now  to  the  term 


H   (we 

n 


■law 


-law 


w 


«W2(1-W)2 


'  X    H   ft)e-2i°'* 


We  have 


U    ^  t   -  w         ^^ 


H   (t)e 

n 

t2ft-i: 


-2iat 


dt    -    (1-w) 


Bjl)e 


-2ia 


t2(t-l)2| 


dt 


t-w, 


00  ,         . 

(1-w)     /    -^ T    dt    +     (1-w)^ 

J      t2(t-l)2 


J    (l)dt 
n 


t2(t-l)2(t-w) 


+   (1-w)' 


,      i  J  (t)-  J  (1) 

^t-lv2     n  n 

^    t    '        (t-1)  (t-w) 


dt 


hG 


where  J  ft)  =  [H  (t)e'^^^^-H  (l  )e"^^'^}/(t-l ).   The  substitution  t  =  (l-x^)'^ 


shows  thsx  for  0  <  w  <  1, 


dt 


t2rt-i)2(t-w; 


W2(l-w)2 


T  tan   r^) 


wsfl-wls 


|lT  -  (1-W)2(1+^(1-W)] 


+  0{(l-w)^) 


Hence 


H  (t  e 
n 


-2iat 


-T J-   (it  -  (1-w) 

t2(t-l)2 


J  t 
n^  ' 


T  dt 


t2(t-l)2 


n  j^(t)-j^(i: 

^^   t2(t-l)' 


2  I'    ~^'^'"'  "'^^  '  dt  +  (J  fl)(l-w)^-  H  (De'^^^'fl-w)) 


+  (1-w)   /   -T TT^ 


n  ■ 


w' 


^n  (l-w)"2-  [1  +^(l-w)} 


+  o((i-w)5/^; 


Therefore 


„  ,  V  -law 

H  (w)e  law 

n^  e 


i  H  (t)e 


-2iat 


la 


w 


i     1_ 

rtw2(l-w)2 


i~)'     \_,         dt  =  -^  ^1  H-  (i  -  ia)(l-w) 

Jt(l-W)2 


^.3   1.    1  2w,   v2l 


H  (t)e 


-2iat 


^1    t2(t-l) 


P  J  (t)dt 
X-  dt  -  (1-w)  /  — T T 

i         t2(t-l)2 


2  P  J  (t)-J  (1)  -2ia  ._   .   „,  ,,w,  .2.   K 


2ia, 


1 


ti(t-i)3/2 


+  2H^(l)e"^^°'  (l-w)-2J^(l)(l-w)^+|Hjl)e-''^^(l-wf 


^7  - 


This  formula  can  be  somewhat  simplified  as  follows: 


-r^^ T  dt 

t2(t-l)2 


H  (t )e      -  H  (l)e 
n n   

ti(t-i)3/^ 


dt 


=  2 


t-1 


H  t  e 
n 


•2iat^ ' 


H  (l)e 
n 


-2ia- 


t2 


2t 


3/2 


dt 


-3/2 
on  integrating  by  parts,  the  factor  (t-1)     being  integrated.   The  prime 

indicates  a  derivation  with  respect  to  t.   Since 


00 

r  (t-i)"*t"3/2dt 


the  last  term  of  the  integral  is  2H  (l)e 


-2ia 


Proceeding  in  the   same  way 


with  the  other  integral  containing  J   (t )  we   find 


H   (w)e 
n 


■law 


w 


law 
e 

I 1 

7rw2(l-w)2 


(^)= 


i  H   (t)e' 


■2iat 


t    -  w 


dt 


la 


(C2) 


3t(l-w)2 


r  (l  +  (|-ia)(l-w)+(|-|ia-|a-)(l-w)^ 


H   (t)e 
n 


-2iat 


L  1 


t2(t-l)2 


dt 


2  1-w 


^+  H    (t)e 

dt        In 


-2iat1 


(t-1)^ 


t2 


.\,.-.,\ 


dt        in 


■2iat 


;t-i)' 


+  0[(l-w)' 
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With  regard  to  d(-w)  the  series  term  may  be  dealt  with  as  in  (C2).   As 
for  the  first  term 


(C3)  -f 


2e 


law 


jrw2(l-w)2    jt(l-w)2 


^^{l+(J-ia)(l-w)+(^-ka-  ia2)(l-wr}  +  0  {  (l-w)^/^} 


Using  (CI)  and  (C2)  we  have 

TTTi       "  i 

2ae  r     (1-it)^      -at 


(C4)      A^ 


t2 


«>  ^    , ,  V    -2iat 
e      "dt   +     > r~~ — T"  dt   , 

Z_j    «  J 

n=l  1 


t2  t-l   - 


2a        h' 

e 
jr 


ni 


(J-ia)i    /    (1-it) 


1  -at 

2  e 


—  dt   + 


t' 


K3_(t. 
t^ 


e        dt 


2n2a2 


:c5: 


la 


n=l 


■la 


H     t   e 
rv 1 


-2iat 


t2(t-l)2 


dt    -   2 


dt 


H   (t)e 
n 


-2iat 


t-iy 


t2 


/I      .     \,      ^  2cte 


2rt^a2 ' 


2e^^      r         dt 


H    (t)e' 


■2iat 


rt        J 


n=l 


i      (t-1 


t2 


:c6: 


A,^    = 


»         ,    ,      1    -,  -at 
A     1-         12,,         A      -u        2ae'       Jl|^2^.    tt      /- (K^(t )  +  ^ile 
(q+21c^  +  2^  )A^+  (--ia)X2-  — ;^ ^3^  a    +  ijt  a  +  J      ^ dt 


2ae  4  2   2,    .    22 


Bit 


la 


dt 


H    (t)e 
n^ 


-2iat 1 " 


n=l 


(t-i: 


t2 


-  k9 


Similarly 

-2iat 


and 


(CIO] 


la 


?<=  \    '     e  /        n ' 


ia       p    E    (t)e 


(C7)  v..      =    ~r     +     7       ~T"  i ~  ^*   ' 

n=l  1 


,  .  -  \T?   f+\    -2iatl 

2e  /  dt        In 


(C8)  ti        =      (2-ic^)^^i   ■    /    ,        n         ;  -   1 

2  2  1  ^  n  J  (^_i)2      I  +2 

n=l  1 


1 


/    .N  ,1     1.        1  2.         ,    ,1      .     s  \  '   ite 

(C9)  ^3     =      (-8+2^^  +  ^    )^^1  +   (i-l<^)^^2       Z_J  ~3^ 


00  Tip   /+  ^    -2iatl 

,    la       p       ,,  E   (tje 

he /       dt       j    n 


n=l 


^1   (t-l)2    I  t2 


1  5 

TJti       00  -f  .^        0°  o    ^-1-  -1   -2iat 

o      ^           r  /I    •4.\2          ,             ^^ — I  ICG      p    S    (tje 

2ae             /  (1-it )        -o-t^,     ,    \    '  e           /     _n__ -,, 

1                        ^               J  t2                                            ^  "         ^            t2(t-l)2 

o  n=l              1             ^ 


2 


;cilj  v^     =     (^-ia)v,    +  —  e"       <    /     -^-i~  e'^-^dt   -   2rt^a2' 


t2 

2iat'| 


^  gH      r        dt  Pn^*)^ 

n.l'      ""       {      (^-1)*     1~~? 


(C12)  v^ 


K-    (t)  +  |i       „^     I  V-1   K   ia       n  Sjtje 


^/-5^=-^4^1>1^i^ 


r^    ,,  ^   -2iat)" 

n' 


jr^'    '-y  'Z..~  J  T^^\~^ 
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Appendix  D 

The  asymptotic  development  of  the  coefficients  \^,\  ,. .. ,v     is  a 

straightforward  matter  i;-^  principle  but  the  detail  is  rather  tedious.   Both 
A^  and  v   can  be  dealt  with  at  the  same  time  by  considering 

1   .     00  I 

O 

with  m=l  and  m=5  in  the  two  cases.   Now 

2ae^     r    (1-it)^  -at,,   2a  l''^    f    e'^M,   1  .,   1  .1   ,  .  P 
o  o       I 


(Dl)  ^^gm^^-XM^ 


+  Trm  (4n  -l)(3m  -2)t^>  dt+0(a^) 


2a-  e         )  -I      ira       Bm    /I        t  ^  ,  5im   /I        n\/l        o\l        ^^      2, 


For  the   series   in  the  coefficients  we  can  use  the  result    (Jones    [17]), 
valid  for  ay  »  1, 

(D2)    H   (y)  ~i   {-^)\   (y)-  -i--  i  ( e^^^^-  i  ) h        +  ^^   (e^i^^^  i 
n^-'  '        jt    ^1+y       n  ^    ]  n^o        12c3ay    ^ 


where 


[m. 


h  =   i     fh     (0)     +   -^    (J    h     (0)     -    h'  (0)1 

n,o       8   [  n^  16a    '•2     n^  n^      ' 


From  this  we  see  that 
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la   p  H  ( t  j  e 

5    /  _n_ -,, 

—  /  r T—  dt 


t2(t-l)^ 


le 


la 


-2iat   h  (t)    h 
:t-l)2   l(l+t)2    ^* 


9h  0 
n 


?  ? 
128a  t 


dt 


la 
e 

2 

rt  a 


dt 


t(t-l)^ 


h    .^°^' 
n,o   128at 


h 


n^ 


+  o(-^) 


a~ 


"jti-ia   00 


-2at 


i_ 

t2 


h  (1-it) 

n 


h 


n,o 


(2-it 


a(l-it) 


9h  (0) 
n 


128a  (1-it 


la   I       9h  (0)       h 
jta   I  n,o    256a  (     ^3' 


on  deforming  the  conto-ur  of  the  first  integral  into  the  line  joining  1  to  1  -  ioo 
and  using 


dt 


t  t-1 


dt 


2     i 


i   2- 


Therefore 
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ia      p    H   (t)e 
e  /       n^ 


-2iat 


"jti-ia       00 


T-  dt   = 


t2(t-l)^ 


-2at 

t2 


h   (1)        .+11 
—    +  -r<Th  (1, 


i_ 
22 


Vn^ 


-h^(i: 


(Dii) 


4<|h:(i)-V(i)+^h(i)l  .V  (!-*-": 


1 


n^    ^     2  n^  Id     n^ 


a 


9h^(0)' 


128a^ 


dt   - 


la 


:rta 


n,o         256a 


h 


Tjti-ia 
e 

2A2 


1, 


^^(i)-i^tVi)-V^)^^2rn 


^"(1) 


Jh'  (1)   +  ^  h   (1)>     + 
2    n^  16     n^      ' 


22h 


^LiO    ,,    _,  i_> 


(1  +   ,     - 


9.22h    (0)" 
128a^ 


jta      1    n,o         256a 


+   0 


n> 


a 


To  evaluate  the  various  terms   in   (Tik)  we  need 


1 

1  -la  -  Tjtl    r 

2  ^ 
■n  e 


a2(l  +  v: 


1   +  2a^l+v     2^    "^  ^2     )8     2(l+v) 


1+v 


.1 
+  0(a   2j 


Hence 
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.      1    .  7 

r2p  "^         ..  15    .         .     "2, 


(D5)  h^d)      -     ^^-^T (l+-^)+0(a      ), 

^  2a2  32a~ 


1.1.  I 

2        -ICC-  ,    JTl  ^  .  ]  J-  Q 

(D6)  h    (0)      =     ^e  ^  (1   +S-  2^    ^°("      )    ' 

^  a^  ^        32a 


1  .      1    . 

2  -ia--r:tl 

:d7)  h^d)     =    -fre        ^     (1*^) 

ija-2 


1  .      1    . 

2  -ia--j^Tti  -2^ 


(D8)  hJO)=  -  ^  e  ^        +  0(a      )    , 

-^  a2 


1  .    1  .  3 

2  -la-TJti  -  ^ 
(D9)               h"(l)      =     -^     e          ^     +  0(a  ^)    , 

•      1    ■  R 

i  -la-Trti  _2 

(DIO)  h  =     i^i^^-T (1  +  iJ)    +  0(a'2) 

1,0  0^2  3^a 


In  our  earlier  paper  it   is   shown  that,    for  n  >  2, 


1 


h  ,  -2ia-^rti  9h     ,(0)      ,  1  /h     / 

(BID      v^)  =  ^^(^-2^)-^— |-Vi^^)-^^^i-2")^°h')^ 

l6TC^a  \  a^ 


-2ia-T7ni 


2       /  32n^a  \a    / 


(-)    ^;a)  =  ^(i-)-"-rvi'^'-iS^'S-^'^°n 


-  5^  - 


(B.3)      VO)      =      -iW^--7|-Vxa)^^^-f-t^ 


(DU)      hJO)      =  Q^ 


Srt^a' 


TTti-2ia  , 


(D15)     h  =     -  -i-  h     ^        +S —  h        ri)   +  0(-^)    . 

n,o  l6a     n-l,o  ^       n-1^  ^      2    ■^ 


-  ?  ^ 


Since  e, (v)   =  2  we  have 
(Dl6)e^(l)   =  2,      e^(0)   =  2,      e|(l)   =  0,      e^(0)   =  0,      e![(l)   =  0,      e^^^  =  i(l  +  j^) 


Formiilae    (Dll)   -    (D15)   are  also  valid  with  h     replaced  by  e   . 


The  corresponding  results   for   s     are 


2      -ICC-  1  rti  _ .         .  __  -  _ 

(D17)  s    ri)      =     ^e  ^      (1-^-i^)    +0(a2)    , 

^  2a2  "-        32a^ 


1.1.  5 


in' 


(D18)  s^(0)      =     ^  e  Ml-  g)   +  0(a  ^)    , 

1 


,2        -la-Trti 


(D19)  s;(l)      =     -  ^e  Ml   -  S)    +  O^'^      )    ' 
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(D20)      s^(0) 


1.1.  i 

^  e  ^     +  0(a  2)    , 


(D21)      s^(i: 


1         .      1    . 
_2      -ia.-r«i 


3 


-^e  ^^^   +  0(a2)    ^ 

1^2 


(B22)      s^^^ 


— I      S 
8a2 


■ia-T3Ti 


(1  - 


lt-)^°(»"'' 


The  formiaae    (Dll)   -    (DI5)  hold  with  s      for  h   . 


■  It   follows   from   {Q,h) ,    (Dl),    (Dl+)-(D15)  that 


(D23)      ^        ^ 

\  =  — I     e        (1 


/  -2ia-7Jti 

-2ia  rr  4 


32a  128a- 


"l6J^3/2 


"a" 


+   0(a'3) 


Prom   (C7),    (Dif),    (DI6)   and    (Dll)-(D15)  we  have 

1 
(D2U) 


^1  = 


2e 


fl   -  — 
n       ^-^        8a 


-Jti-ia 


^).s- 


128a 


jt  a 


-2ia-Trti 
.+3^._lli-+e 

^^     512a2  3 

I6rt2a 


+  0(a'3) 


From   (CIO),    (Dl),    (D17)-(D22)   results 


i     1    . 
2a2      ^"1 


(D25)      V      =  — r  e 


i-a,JtL^  891  ,^ 


-2ia 


ifjta 


-51^^ 


-2ia-Tni 


I6rt2a 


+  0(a' 


56 


Proceed  in  a   similar  maniier   for  the  coefficient    Xp^Hp  and  Vp.      We  have 


J\-  e"°-*dt   -   2A^    =|^^^-|mi-  ^(|n-l)t  +  -^m(im-l)(|n-2)t2[>dt 


(D26; 


2jt^a2    +   0(a   ^) 


Also 

^       J      (t-l)2     I  t2  J  ^  {  t2(t-l)2 


and 


H^(t)[   ^-2iat  g.^ia    Pel^^^  J_^>L     _V^     SA  .. 

:t-l)2  /      -^     (t-l)2      (1+t)-^  ^       at^ 

1  ^  2(l+t)^  ^ 


.2    J    t(t-iFp^°     '^^'1    ^'-    ^    t2(t-i)i     ^^2) 
1  ^  -^  1 
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-jri-ia 


2e 


2  t 
ir  a. 


{h,(i: 


k(i)}- ^„- (kd)- <(i)  +  't>':(i)) 


ij.  n^      -'      32cx   '■li  n^ 


n  n 


(D28; 


22a     ^'° 


, .    ia  9h    (0)        ih^   _  h 

'  a 


in  the  same  way  as  was  employed  in  deriving  (d4).   Use  (D4)  for  the  first 
term  in  (D27)  and  then 


2e 


la 


(t-l)2 


dt      \^^^' 


-2iatl t 


t2 


1    1  .  . 

2ia'^   4 
2 


^(i)-tK(i)-k(i)) 


Vn^ 


(D29) 


+  ^  {h"(l)-  ih'(l)+T%  h  (1)1+ (1-  f-, 

2  '■  n^    2  n     16  n^   -"    a    ^    i|a 


32a 


2%^o 
a 


i_ 

9.2^h  (0)' 
n^ 


128(x' 


jta      ^2' 
a 


The  combination  of  (C5),  (D26),  (D29),  (D5)-(D15)  gives 

(D30)  I   1  • 

1   .  s,    i*x    if"^^,,  ,  i     1  ,  _i_,  ,  i  -2ia,,  ,  i   e 


-2ia-Tjri 


^2  "  (2"^°"^'  — ~  ^  ^'^'^ 


^      32a2   ^20^ 


3^";^"    (1"I6^^- 


.1 
le^t^a 


•)  +  0(a 


-2. 


From    (C8),    (D29),    (DI6),    (D11)-(D15)    follows 


1      1 


,1 


•    ICC      ),  .^■2     TJti-ia  .  ,  ^ 


-2ia--rjti 
4 


'°31)«.  =  (5-1.^-1;;:^=!-  e'  (^■.^^^,^3  3 


+  0(a"^). 


l6n' 


a 
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In  addition,  (Cll),  (D26),  (D29),  (D17)-(D22)  give 

(D32)  I  -\t,±  -2ia-^7ti 


/1.x    ki  e ,^      5i   15    Hi  X   ie ,^      kl±  ,   e .  ^, 


T  ^  2  2 
16«  a 


The  coefficients  X^.m--,  and  v^  can  "be  determined  in  a  similar  manner. 
3  J      J 


Thus 


3                           °°        ,    ^     1  3 

1,     i  Q            i  i        /-■    K  (t )  +  rmi  ,  i    i    o            i    i 

^     2   2  ,    .      2    2         /        111            2  -at , ,  4  p   2             p-  p 

-  j:  a    +  imjt  a''  +    / e        dt    =  -jr  a    +  imir  a'^ 


°  t^ 


-at        ,      ,  .      ,  .  -|^ 


(D33)  +    /— -  {-^(^-l)+Y|(im-l)(|m-2)}dt   +  0(a 


t2 


1   ^  -5 

k  2   2  f^  ^  Sim        Sin    /I      -,  \,      im    ,1     .,  . /■!     ^o         ^,      2- 
=  :^rt  a     |l+-j^-  --^(^-1)+  r(^-l)(^-2))    +  0(a      . 

■^  l6a  32a 


Also 

00  f^T    /^  \    -2iati''  o  °°TT    /j- ^    -2iat 

ite^      r   __dL  IVllL I  _       _l6a!  ^ia    f^n^^^"  ^^ 

3^      {      (t-l)*!  t*  J  3^     ^      {        t^(t-l)* 
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The  first   integral  is  given  iDy   (D4)   and  the   second  by   (D28).      As   for  the  third 


he 


xa. 


°   -2iat      H   (t) 

^ rW 

(t-l)2  t^ 


dt 


U> 


la 


-2iat        h    ft) 
e  J    n-     ' 


3« 


;t-i)2    l(i+t)2 


h^(t)_^jyo 


A  1 

(1+t)^      4(l+t)^ 


dt 


l6ai 


la 


dt 


3rt 


t(t-l. 


n,o   ^  128at    f  ^        ,^2    J      t"(t-l)2 


3rt 


2e 


Tjri-ia 


K'l'-iv^ 


l6a6 


3« 


ia   f  9h^(0)        ih^^^ 


^,0"^     256a     "^     2a 


h 

+  o(-a-) 
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Consequently 


lie 


la 


dt 


H   (t)e 


-2iatl 
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3« 


:t-i)^ 
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8a2      4"i-i" 

,   2 


h    (l)+f- 
n^  iicc 
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h    fl) 
n  ■ 
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(D311) 


32a 


>;(!) 


k(^' 


+  -^  h  (1 

16     n 


—  h        (1 
a      n,o 


31  ■ 
ija' 


128: 


i  \(o: 


h 

+  0( — , 
^  a 


By  means   of   (C6),    (D33),    (D3^),    (D5)-(D15)  we   find 


,1     1.        1  2. 
(o  +-ia  +  Ta 


2     Trti 


^3   "    ^8  "2"^  '2^    ^^  •    ^2 


,1      .    ^,        8x"     V     ,,^  3i  ^     3     ^    3i      s 
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^  ■  32a2  •  ^28a3 ' 
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-)    +  0(a      )    . 
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A  consequence   of    (C9),    (D3^),    (Dl6),    (Dll)-(D15)    is 


^     1    .     .  -2ia-7ni 


(D36)  ^,    =    (i,iia.^2^^^    ^    ^i  -iaK-^e^"""(l-Tlr--^^^     ^  Wo    ) 


3/ 

From    (C12),    (D33),    i'Q'ih) ,    (D17)-(D22)   it   follows   that 

,1,1.     ^1   2.           ,1      .     ,           8a^     1^"^,,      151        k'^  l-^i    , 

V      =    (Q+-ia+^    )v     +(--ia)v    -  — r  e        (1+^—  '-  +  —^~] 

^^^  3^2  ht      32a^      128a3 


(D37) 


-2ia--rjri 
hp.     -2ia,^      ';J3i  ,  e x        ^,    -1. 
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